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Most of the analysis of the Cosmic Microwave Background relies on the assumption of 
statistical isotropy. However, given some recent evidence pointing against isotropy, as 
for instance the observed alignment of different multipoles on large scales, it is worth 
testing this assumption against the increasing amount of available data. As a pivot model, 
we assume that the spectrum of the primordial perturbations depends also on their 
directionality (rather than just on the magnitude of their momentum, as in the standard 
case). We explicitly compute the correlation matrix for the temperature anisotropies in 
the simpler case in which there is a residual isotropy between two spatial directions. As 
a concrete example, we consider a different initial expansion rate along one direction, 
and the following isotropization which takes place during inflation. Depending on the 
amount of inflation, this can lead to broken statistical isotropy on the largest observable 
scales. 
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Qh| 1. Introduction and discussion 

O I The release of the WMAP 3-year results^ has now placed our knowledge of the 

' largest scales of the CMB on a firm footing. An intriguing aspect of the maps 

. produced by WMAP is the presence of a number of seemingly related anomalies 

on the largest scales. They comprise the lack of power in the lowest multipoles, 
L J the alignment of the power in a number of the lowest multipoles and an apparent 

, asymmetry in the maps. The lack of power and its possible origin has been addressed 

' copiously in the literature. ^'^'^ The difficulty of calculating the posterior probability 

of the measured values of the quadrupole and octupole has led to a spread in the 
claims of the significance of the effect which now appears to be marginal.'^ 

The significance of the anomalies which represent a breaking of the statistical 
isotropy of the maps appears much stronger. ^'^ Indeed, if the signal is statistically 
anisotropic then much of the information content of the maps lies in the anisotropic 
correlation of the spherical harmonic modes with (af„iOf/,„/) 9^ See'Smm'- Obvious 
explanations for the observed anomalies are that an unknown systematic is not 
being accounted for in the analysis or that foreground signals, which are naturally 
anisotropic in nature, are not being modelled and subtracted properly. However a 
conclusive explanation along these lines has not been put forward yet. 



*Talk given by Carlo R. Contaldi at the Eleventh Marcel Grossmann Meeting on General Rela- 
tivity, Berlin, July 2006. 
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Another, albeit more speculative, possibility is that the explanation lies in more 
fundamental aspects of the cosmological model. It seems advantageous therefore to 
attempt to model the possible anisotropics in the sky from first principles from a 
primordial perspective with a view to constrain departures from the standard cos- 
mological picture. Here we assume that the broken statistical isotropy is imprinted 
in the primordial spectrum P$ (k) of the cosmological perturbations, leading to a 
specific and falsifiable pattern for the covariance matrix of the CMB fluctuations 0. 

The present discussion is organized as follows. We first compute the general cor- 
relations that are obtained in this case. We then present a simple model which gives 
rise to the above correlations. The background evolution of this model is axisymmet- 
ric, with an expansion rate which is initially much faster along a single direction. 
For simplicity, we assume isotropy in the remaining two directions. A period of 
inflation than leads to full isotropy. Such a situation may emerge in models with 
extra dimensions (as for instance in string theory). In this case, one faces the diffi- 
cult problem to explain why only three directions undergo cosmological expansion, 
while the other remain stabilized on a small size.^ It is plausible to imagine that 
the expansion of the three (now) large dimensions did not start precisely with the 
same rate, but that isotropy was only reached afterwards in the first few e-foldings 
of inflation. With this approach, matter and radiation perturbations propagate in 
an isotropically expanding background. The breaking of the statistical isotropy is 
therefore entirely due to the initial conditions of the model (as opposed to a full 
anisotropic evolution of perturbations^'^*^). We conclude by showing the spectrum 
of perturbations of scalar modes having the momentum in the anisotropic direction 
(as we discuss, the computation is considerably simpler in this case). 

Details of the computations, and the comparison with the observations, will be 
presented elsewhere. 

2. CMB covariance matrix 

We want to calculate 

C«'mm' = {ai^a},^,) = j dnpdnp'{ST{p,rio,Xo)6T{p\r]o,xo))Yg*^{p)Yg,^,{p'), 

(1) 

where ST is the temperature perturbation in the direction p, as measured by an 
observer at position xq and at time 770 ; it can be expressed as 

^<i>(k,r?OA(A:,k.p,r;o)e'''''°, (2) 

We assume that the anisotropy is only imprinted in the primordial curvature 
power spectrum. Specifically, {2Tr/k)^S^(k — k')P^{k) = ($(fc)$*(/c)) depends also 



''This is in contrast to the opposite approach of constraining general modifications to the covari- 
ancel'^ such as (a^,„a^/„/> oc <5«/<5„„/(C^ + m^). 
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on the orientation of the wavevectors, and not only on the magnitude. In con- 
trast, the calculation of the later time transfer function Af(fc,77o) for the radiation 
perturbation proceeds as the standard case by the line of sight integration of the 
Einstein-Boltzmann system. Note that the transfer functions are isotropic in that 
they only depend on the magnitude of the wavevector k. However the final integra- 
tion over wavevectors k must be modified due to the directionality dependence of 

For definiteness, we assume a residual symmetry along two directions. We denote 
by ^ the cosine of the angle between the wavevector and the axis of symmetry. 
Working the expression ([T]) out, we find 



y {£ + m)\{£' + m)\ 

X J ^A,{k,m)A,{k,rjo) f^dCPnOPl'iOP-fikA) (3) 

where A; are the coefficients of the decomposition of the transfer function into 
Legendre polynomials). 

In the isotropic case, P.^ (fc, ^) = F$ (fc) , we recover the standard result 

Cu-tm- = -5u'5„^,n' j y {k) I A, (fc, ryo) ? (4) 

In general, we obtain a correlation between different multiples. Here, we focus on a 
axisymmetric model, for which P^(k,—^) = In this case, the correlation 

is non- vanishing whenever the difference between I and I' is even (as it follows from 
the parity properties of of the associated Legendre polynomials). 



3. Background evolution 

As a specific example, we consider a homogeneous but anisotropic background, with 
an expansion rate along the x-direction different from the other two, 

ds^ = -dt^ + a {tf dx^ + b {tf {dy^ + dz^) . (5) 

A scalar field in this background leads to the equations 

ij + 3ij2 = y, 

iH^ -h^ = ^4,^ + V, (6) 

where H = {Ha + 2Hi,) /3 and h = {Ha — Ht) /VS (we work in units of Mp — 1, 
while prime denotes differentiation with respect to 4>). The first two equations are 
identical to the ones obtained in the isotropic case, in terms of the "average" expan- 
sion rate H . The third equation, appearing as a "modified" Friedmann equation, 
can then be used as an algebraic equation for the difference h between the two 
expansion rates. 
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Fig. 1. Evolution of the Hubble rates and the difference in e- foldings during the anisotropic 
and the following isotropic inflationary stages. We have chosen a chaotic inflationary model, with 
quadratic potential V = m^ip^ /2 with initial conditions 0o = 16 , /lo = 10* ™0 ■ This leads to 
about 60 e-foldings during the isotropic inflationary stage. 



We consider a scalar field potential, and an initial scalar field value capable of 
sustaining inflation, and a large initial anisotropy, /iq 3> Vq S> 0o ■ The system 
then undergoes a quick anisotropic inflationary stage (with both a, 6 > 0), during 
which (f> is practically static, i.e. V Vq, and the two Hubble parameters behave 
as Ha — l/t and Hb ~ 0. To follow the evolution of h, it is useful to combine 
eqs. ([6]) into h + 3H h = 0. This leads to a very rapid decrease of h during the 
anisotropic stage, as the average rate approaches the standard (slow roll) value 
H ~ . The timescale for the decrease is ^iso = \/3/Vb (we can also see this 

by requiring the continuity of Ha between the early l/t behaviour and the later 
constant value). From this moment on, h can be neglected, and we have a second 
stage of standard (isotropic) inflation. Figure [1] shows, for a specific choice of initial 
conditions, the evolution of the Hubble parameters along with the difference in the 
number of e-foldings between the expansions in the x and y, and z axes, defined as 
A A'' = J* dt {Ha — Hb). This increases during the first stage and remains constant 
after tiso • 

4. Perturbations 

The standard FRW background is symmetric under spatial transformations. It is 
therefore convenient to classify the perturbations according to how they transform 
under spatial change of coordinates: modes which transform differently are not 
coupled at the linearized level. One then finds that 3 physical modes are present. 
They are contained in one scalar perturbation, and in the two polarizations of a 
tensor mode, which are also decoupled from each other. Also in the present context 
there are three physical modes. Two of them are coupled to each other, while the 
third one is decoupled, due to the residual symmetry in the y — z plane. The first 
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Fig. 2. Power spectrum for longitudinal scalar perturbations. The initial conditions for the back- 
ground are as in figure [T] while </> = 3 at the moment shown here. 

two modes decouple when the universe become isotropic, and they become a scalar 
and a tensor polarization. 

Therefore, we expect that the coupling is proportional to the difference between 
the two Hubble rates. We also expect that the coupling vanishes when the momen- 
tum of the modes is along the privileged direction x . We denote these modes as 
longitudinal. Indeed, the directionality of a mode always breaks the full 3d isotropy; 
when the momentum is oriented along the anisotropic direction x we have the same 
symmetries as in the standard case. Therefore, the computation is much simpler for 
these modes. Here, we focus only on this simpler case. The complete computation 
is given in.^^ 

We proceed as in the standard case, by writing down the most general scalar 
perturbations of the metric and of the inflaton field, and by choosing a gauge which 
completely fixes the freedom of coordinates transformations (see^^ for a review). 
In the longitudinal case, the linearized Einstein equations for the perturbations 
can be reduced to a unique equation in terms of a single variable Q , which is the 
generalization of the Mukhanov-Sasaki variable, and which obeys the equation 



where k is the comoving momentum, oriented in the x direction. This equation 
coincides with the standard one when Ha = Hi, = H. We note that the scale 
factor a is entering in the physical momentum k/a, since the mode is moving along 
the x-direction. The initial conditions are given as in the standard case, since k/a 
dominates the last term of eq. ([7]) at early times, leading to a standard adiabatic 
initial vacuum. 



Q + iHa + 2Hi,)Q+ —+M^ g = 
M = V" + 2<j)V'/Hb - 0V {2H^) + 2Ha<j)^/Hb + <j)' 



i2 



(7) 
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In figure 21 we show the power spectrum towards the end of inflation, kiso is 
the conioving momentum of the modes which left the horizon when the universe 
became isotropic. Modes with k > kiso leave the horizon during the later isotropic 
stage of inflation. The standard result is recovered for this modes, with the typical 
spectral index of massive chaotic inflation. However, modes with smaller momenta 
left the horizon when the universe was still anisotropic, and a nonstandard result is 
found in that case. 

Acknowledgements 

We thank Joao Magueijo for useful discussions. The work of A. E. G. was supported 
by the Hoff Lu Fellowship in Physics at the University of Minnesota. The work of 
M. P. was supported in part by the Department of Energy under contract DE-FG02- 
94ER40823, and by a grant from the Office of the Dean of the Graduate School of 
the University of Minnesota. 

References 

1. Spergel, D. N., et al. 2006, 'arXiv:astro-ph/0603449' 

2. Bielewicz, P., Gorski, K. M., & Banday, A. J. 2004, Mon. Not. Roy. Astron. Soc, 355, 
1283 

3. Bridle, S. L., Lewis, A. M., Weller, J., & Efstathiou, G. 2003, Mon. Not. Roy. Astron. 
Soc, 342, L72; Schwarz, D. J., Starkman, G. D., Huterer, D., & Copi, C. J. 2004, 
Physical Review Letters, 93, 221301; Slosar, A., Seljak, U., & Makarov, A. 2004, Phys. 
Rev. D, 69, 123003. 

4. Contaldi, C. R., Peloso, M., Kofman, L., & Linde, A. 2003, JCAP 0307, 002 (2003). 

5. Efstathiou, G. 2004, Mon. Not. Roy. Astron. Soc, 348, 885; de Oliveira-Costa, A., 
Tegmark, M., Zaldarriaga, M., & Hamilton, A. 2004, Phys. Rev. D, 69, 063516. 

6. Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2006, 
arXiv:astro-ph /0605135 Eriksen, H. K., Hansen, F. K., Banday, A. J., Gorski, K. M., 
& Lilje, P. B. 2004, Astrophys. Journal, 609, 1198; Eriksen, H. K., Banday, A. J., 
Gorski, K. M., & Lilje, P. B. 2005, Astrophys. Journal, 622, 58; Hansen, F. K., Banday, 
A. J., & Gorski, K. M. 2004, Mon. Not. Roy. Astron. Soc, 354, 641; Hansen, F. K., 
Cabella, P., Marinucci, D., & Vittorio, N. 2004, ApJL, 607, L67; Jaffe, T. R., Banday, 
A. J., Eriksen, H. K., Gorski, K. M., & Hansen, F. K. 2005, ApJL, 629, LI; Land, K., 
& Magueijo, J. 2005, Physical Review Letters, 95, 071301; Land, K., & Magueijo, J. 
2005, Mon. Not. Roy. Astron. Soc, 357, 994; Land, K., & Magueijo, J. 2006, Mon. Not. 
Roy. Astron. Soc, 367, 1714; Ralston, J. P., & Jain, P. 2004, International Journal of 
Modern Physics D, 13, 1857; Vielva, P., Martinez-Gonzalez, E., Barreiro, R. B., Sanz, 
J. L., & Cayon, L. 2004, Astrophys. Journ al, 609, 22. 

7. Magueijo, J., & Sorkin, R. D. 2006, arXiv:astro-ph/0604410) 

8. R. H. Brandenberger and C. Vafa, Nucl. Phys. B 316739T (1989). 

9. Bunn, E. F., Ferreira, P. G., & Silk, J. 1996, Physical Re view Letters, 77, 2883 

10. Gordon, C., Hu, W., Huterer, D., & Crawford, T. 2005, |arXiv:astro-ph/050936l| 

11. Contaldi, Giimriikgiioglu, Peloso, work in progress. 

12. V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 
(1992). 

13. V. F. Mukhanov, JETP Lett. 41, 493 (1985) [Pisma Zh. Eksp. Teor. Fiz. 41, 402 
(1985)]; M. Sasaki, Prog. Theor. Phys. 76, 1036 (1986). 



